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FOREWORD 


This report entitled "Development of Response Models for the Earth Radiai >n 
Budget Experiment (ERBE) Sensors" consists of the following four parts. 

Part I, NASA CR-178292, is entitled "Dynamic Models and Computer Simulations 
for the ERBE Nonscanner, Scanner and Solar Monitor Sensors". 

Part II, NASA CR-178293, is entitled "Analysis of the ERBE Integrating Sphere 
Ground Calibration". 

Part III, NASA CR-178294, is entitled "ERBE Scanner Measurement Accuracy 
Analysis Due to Reduced Housekeeping Data". 

This is Part IV, NASA CR-178295, entitled "Preliminary Nonscanner Models and 
Count Conversion Algorithms". 
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1. INTRODUCTION 


This report entitled "Development of Response Models for the Earth 
Radiation Budget Experiment (ERBE) Sensors" consists of four parts. This 
part, Part IV, NASA CR-178295, is entitled "Preliminary Nonscanner Models 
and Count Conversion Algorithms". 

The remaining parts are as follows. 

Part I, NASA CR-178292, is entitled "Dynamic Models and Computer 
Simulations for the ERBE Nonscanner, Scanner and Solar Monitor Sensors". 

Part II, NASA CR178293, is entitled "Analysis of the ERBE Integrating 
Sphere Ground Calibration". 

Part III, NASA CR-178294, is entitled "ERBE Scanner Measurement 
Accuracy Analysis Due to Reduced Housekeeping Data". 

This document discusses in greater detail the preliminary nonscanner 
models introduced in Part I of this report. The mathematical modeling 
equations and count conversion algorithms are discussed in greater depth. 
The purpose of this document is therefore to define two count conversion 
algorithms and the associated dynamic sensor model for the M/WFOV non- 
scanner radiometers. The model is required to provide and update the 
constants necessary for the conversion algorithms, though the frequency 
with which these updates needed to occur was uncertain. Since the count 
conversion equations themselves are functions of the sensor temperatures 
and impinging radiation in which these temperatures are transmitted to- 
gether with the sensor radiometric data, other parameters updated from 
the senosr model were suggested to be as infrequent as once per month after 
instrument performance had been validated. The final selection of an al- 
gorithm and its associated sensor model was dependent on the completion of 
(1) sensor characterization and calibrations and (2) computational simula- 
tion and analyses of various measurements and data reduction processes. 



2. NARRATIVE DESCRIPTION 


2.1 Approach - The following approach was pursued: (1) To develop 

a mathematical model of the sensors which accounts for the conversion of 
2 

irradiance (W/m ) at the sensor FOV limiter to data (counts). (2) To 

derive from this model two algorithms for the conversion of data (counts) 

2 

to irradiance (W/m ) at the sensor FOV aperture. (3) To develop measure- 
ment models (which account for a specific source together with the sensor) 
for the 

(a) sensor characterization and calibration, in order to empirically 
determine constants in the sensor model and conversion algorithm; 

(b) pre/post launch calibration, in order to check and update (if 
necessary) constants in the sensor model; and 

(c) inflight measurements, in order to properly apply constants to 
the conversion algorithm. 

Of the two algorithms, one is of the gain/offset type that would be 
preferred because of its computational simplicity provided that it is suf- 
ficiently accurate. The other algorithm is of the Kalman filter type that 
may otherwise be required. The mathematical models of the sensors and 
count conversion algorithms are presented in this document. 


2.2 Sensor Models - The nonscanner sensors are active cavity radio- 
meters with medium and wide f ield-of-view (M/WFOV) limiters and shortwave 
(SW) and total (T) responses (see, e.g. Fig. 1). The SW response is ob- 
tained by the inclusion of a fused silica dome filter over the cavity 
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aperture. The dynamic model for the T sensors has 8 themal nodes and for 
the SW sensors has one additional node for the filter. In addition, the 
model accounts for the active electrical heater (H) and three temperature 
probes (P) as well as for the signal amplification and analog-to-digital 
conversion. 

Irradiance inputs, which are treated as another thermal node, are 
partially incident directly on the cavity aperture (node 8), partially ab- 
sorbed and reflected by the base (5) and FOV limiter (6), and for the SW 

i ' .i 

sensor also absorbed, reflected, and transmitted by the filter (9). The 
absorbed irradiance raises the temperature of nodes 5, 6, and 9 and causes 
unwanted thermal radiation toward the cavify (1) and heat conduction through 
nodes 4, 5, 6, 8 and 9 to the cavity (1) and heat sink (2). Some of the 
reflected irradiance from nodes 5, 6 and 9 may scatter into the cavity (1). 

The unwanted as well as source irradiance that is absorbed by the cavity 
walls increases its temperature. The electrical resistance temperature 
probe wound around the cavity walls senses this temperature change and 
activates an electrical feedback circuit which, in turn, supplies power to 
an electrical resistance heater that is also wound around the cavify walls, 
(see Fig. 2). The feedback electronics (see Figs. 5 and 6) are designed to 
maintain a constant temperature difference (approximately 1°C) between the 
cavity and the reference cavity. The heat sink, and hence the reference 
cavity, is temperature controlled to approximately 40°C. The electrical 
voltage required to maintain this constant temperature difference is (ideally) 
proportional to the heat inputs to the cavity. The voltage signal is 
amplified, converted into digital counts, and transmitted together with house- 
keeping data, which includes the temperature of the FOV limiter (6) and the 
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base (5) and copper (2) heat sink. 


2.3 Conversion Algorithms - Two types of conversion algorithms have 
been developed from the sensor model. The steady-state linearization (SSL) 
algorithm, which is of the "gain/offset" form, requires fewer on-line comp- 
utations; however, it may have larger errors due to the sensor's transient 
behavior. The modified Kalman Filtering (MKF) algorithm is more accurate 
under transient conditions, but requires more on-line computations as well 
as filter design development. 

Both algorithms (implicitly or explicitly) account for temperature dif- 
ferences at different nodes, temporal changes in the node temperatures (e.g., 
dome filter, FOV temperatures), reflection and emission from FOV limiter in- 
to the cavity, active cavity heater behavior, electronic feedback control, 
etc. However, the final accuracy of these algorithms (in comparison to 
other algorithms) can only be ascertained by testing and simulation under 
realistic conditions. Thus, ground calibration data must be used for testing 
and development . 

2.4 Measurement Models - Figure 4 presents a functional flow diagram 
of the sensor, its model and count conversion algorithm, and the error 
analysis process that is used to 

(a) evaluate the sensor model (error e^), 

(b) evaluate the count conversion algorithm (error e^)» 

(c) adjust model and algorithm parameters to minimize errors and 

(d) establish confidence intervals for the converted data. 

This analysis must be performed with sensor characterization and calibration 
data. Characterization and calibration sources provide a variety of irradi- 
ance inputs. At the same time, unfortunately, stray heat inputs leak 
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i 

(radiatively and conductively, directly and indirectly) into the sensors 
from the sensor/ source environment. If a sensor is calibrated in a given 
environment but is used in a different environement , errors are inevitably 
introduced unless the effects of the calibration environment are properly 
accounted for. Thus, it is important to account for these stray heat in- 
puts when calibrating the sensor to empirically update constants for the 
model and algorithm. Figure 5 shows a block diagram of the procedure and 
approach that will be used during pre/post calibrations. The nominal 
values of the model and count conversion parameters have been established 
in the ground calibration phase. The approach in this phase is to: 

(a) determine whether the sensor response has changed, 

(b) determine whether the calibration sources have degraded, 

(c) update the sensor and source parameters to adequately reflect 
the new (usually degraded) sensor response, 

(d) update the count conversion parameters and 

(e) update the confidence intervals on the measurements and count 
conversion. 

Thus, current inflight calibration data is compared to the corresponding 
previous calibration data to determine if any change in overall response has 
occurred. If significant differences are found then sensor and source model 
parameters are varied systematically, and the current data is now compared 
to simulated data with updated parameters. When updated parameters are ob- 
tained error analyses are performed to obtain new confidence intervals, as 
well as updated count conversion parameters. 

Figure 6 shows the functional flow of the measurement process including 
the sensor /environment interaction, routine data processing and the inter- 
action of flight calibration processing with routine data processing. 
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3- SIMULATION INPUTS 


The following Input Data Table I outlines the constants and parameters 
which needed To Be Determined (TBD) for the corresponding element numbers 
and locations for use in the nonscanner sensor model (see sections 5*1 and 
5.2). 


Input Data Table II further lists the contants and parameters which 
needed To Be Determined (TBD) for the nonscanner model equations. 


INPUT DATA TABLE I 

(To Be Determined For The Sensor Simulation) 







ELEMENT 

NUMBER 

OR LOCATION 





Constant or 
Parameter 


n 

1 

e 

2 e 

3 f 

3 a 

3 

e 

A 

c 

m 


El 


8 b 

KB 

9 

a 


P S 

TBD 

BSj 

EH 



TBD 

TBD 

TBD 



TBD 



B 

TBD 

TBD 


TBD 

P L 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

S 

p s 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

s 

P L 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

■ 

TBD 

TBD 

i 

TBD 

TBD 

TBD 

TBD 

up 

TBD 

P S 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

4 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

€ s 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 


TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

T 

■ 

■ 

■ 




■ 

■ 


■ 

■ 

■ 

■ 


TBD 

TBD 

T 

■ 

■ 

■ 

i 



■ 

■ 




■ 

■ 


TBD 

TBD 

k 

TBD 

■ 

■ 

TBD 

TBD 

■ 

■ 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

— 

TBD 

TBD 

TBD 


TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

! TBD 

TBD 

TBD 

V 

TBD 



TBD 

TBD 



TBD 

TBD 

TBD 

TBD 


TBD 

— 


TBD 

m 

Y 

TBD 

■ 

■ 

TBD 

TBD 

l . 

■ 

■ 

TBD 


TBD 

TBD 


BE 


TBD 

■ 

C 

TBD 





TBD 

TBD 


■ 

TBD 

TBD 

TBD 

TBD 


TBD 


TBD 

5 


6 






































































































































































































IMPUT DATA TABLE II 


Vi 

L i-k 

Vi 

Vi 


TBD 

k 

m 

1 , 

2, 

• • • 9 

11 

or 12 

TBD 

i 

m 

1 . 

2, 

• • • t 

11 

or 12 

TBD 

k 

m 

1 , 

2, 

• • * 9 

11 

or 12 

TBD 

k 

m 

1 , 

2, 

• • • » 

11 

or 12 

TBD 

k 

MX 

1 , 

2, 

3 




(dependent on Sensor operation schedule) 


Fj_ i - TBD 


k - A, 5, 6, 7, 8, 11 
j & i - 1. 8,, 8 , 11 

(with filter dome) 

or 5, 6, ll b , 12 
- 1. 8 , ll a 

(without filter dome) 

or 5, 6, 8 t . Ilf 12 


F j-i (K-L)° TBD 


■'FOV(t) 


TBD 


(on account of the variety of radiant sources, 
the radiant sources must first be defined before 
determining E F()V(t) .) 


Housekeeping Data (HKD) & Necessary Information 
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4. SIMULATION OUTPUTS 


After incorporating the constants and parameters from Input Data Tables 
I and II, the simulation model provided the following outputs. These out- 
puts became inputs for use in further computations. 

(a) All parameters defined in Section 5 for use in the Data Validation 
Algorithm for data products 

(b) Noise, bias and time constants and some input variables for use in 
the Error Analysis problem 

(c) Irradiant flux values at FOV of sensor for use in the Inversion 
Analysis Algorithm. 

5. PROCESSING METHOD/ ALGORITHM/ EQUATIONS : 

The equations for the radiative and thermal models of the M/WFOV NS 
sensors with and without a filter dome are given in Equations 8 and 9 in 
section 5.2. The coefficients A, D, B^, B and b in these equations are 
computed from Input Data Table I and the variables defined in this section. 

Equations 8 and 9 after having defined coefficients A, D, B^, B and b, 
radiometric count readings and HK data, are used in the Count Conversion 
Algorithms for the analyses of both steady-state and transient responses 
of the sensors. 

A simulation of the sensors will be developed using the parameters de- 
fined in this section. This simulation will be used as shown in Figures 
5 and 6 during ground and flight calibrations to update model parameters, 
count conversion parameters, for error analyses and the selection of 
algorithms . 

The following pages show how the coefficients A, D, B^, B and b of 
Equations 8 and 9 are defined. 
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5.1 NS SKNSOR MODEL 


The constants, parameters and variables which appear in the following 
equations are taken from Input Data Table I or defined by the additional 
equations in Input Data Table II. 

The indices of variables are denoted by the node surface identifier 
shown in the attached figures, e.g., 9a, 3d, etc. 

Exchange factors ; 

E j-i " F j-i + (E P l) r P K F j-i(K-L) (1) 

where i, j = identifiers for each surface 

L »• 0, 1, 2, ... number of specular image surfaces. 

K - 1, 2, 3, ... number of specular real surfaces. 

Variable 1 : 

P j-i “ (<S ji " p j e j-i* 
where i, j = identifiers for each surface 

Variable 2 : 

V ■ s Vj p j-i <3) 

where k, j, i = identifiers for each surface 
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Variable 3 & A: 


With filter dome: 

5 ’ i 1 - *9 \j9 a \»9 b ] 

where j - 1, 8 fe , 8 t> 9 & 
m ■ 5, 6, 9^» 10 

Without filter dome: 

G = f 1 - T 9 H 9a j9 a H 9 b m9 b ] 

where J ■ 1, ^ 9 g 

m = 5, 6, 8 t » 9^» 10 

Variable 5 : 

\ • <w >k 
where k * node number 


Variable 6: 


*ki 


hc-i + 1 + L i-k 

\-l \-i \-i K i \-i 


where k, i ■ node number 
Variable 7 : 


11 


<- 4- K l K 2 K 3 )/M l 


(A) 


(5) 


( 6 ) 


(7) 


10 



12 R lA - 


A 13 ' K l K 2 K 3' A 1 10 " K 3 


all other 1st row elements In the A matrix are zero (see Eq. 8), 


21 


'22 


R 21 M 2 


'( R 21 + R 24 + " R 23 + 4 ) ‘ (0 


■23 R 2J M 2 


24 R 24 M 2 • 


27 R 27 M 2 


all other 2nd row A matrix elements are zero. 


32 R 32 M 3’ 


33 


R 32 M 3 


all other 3rd row A matrix elements are zero. 
1 


* 42 ' R 42 M 4 


44 


' "( R « + R 45 + R 48 ) * l R 4 ) 


A « ' « 45 m 4 ’ 


48 R 48 M 4 


all other 4th row A matrix elements are zero. 

. - - -1- A . JjL. + -L- + _L. + /_L 

55 \ R 54 R 56 R 58 R 51 1\57 


'54 R 54 M 5 • 


)) * (k) 


57 


R 57 M 5 


11 


"56 - R 56 M 5 • *58 ’ R 58 M 5 • *59 - Rj, Mj * "51T R,,^ 

all other 5th row A matrix elements are zero. 


(with filter dome 


65 R 65 M 6 ’ 


A, , = - 


66 R 65 M 6 


all other 6th row A matrix elements are zero. 


72 


R 72»7 ’ * 77 '(' R 72 M 7 ' R 75 M 7 ) 


•75 R 75 M ? 


all other 7th row A matrix elements are zero. 


84 R 84 M 8 ’ 


85 r 85 m 8 


a 88 " -(4 + 4 + * (y (with fiiter a ° me> 


or 


(without filter dome) 

V R 84 R 85 ' W 


8 , 


(with filter dome) 


'll R 118^8 

all other 8th row A matrix elements are zero. 
For a sensor with filter dome only : 


811 R 95 M H 


118 • R M 
118 8 


A n n = (“ r^s + Rn 8 ) (hh) 

all other 9th row A matrix elements are zero. 
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Considering the box beam as an additional node designeated "node 9", 



all other 10th row elements are zero. 


Variable 8: 


The subscripts S and L denote the shortwave and longwave properties 
of surfaces in the enclosure, respectively. That is, the radiative 
properties in the bracket with subscript S are defined by shortwave 
characteristics of a surface. These properties are taken from Input 
Data Table I. 

To avoid the repetition of writing, we write as the following 
example: In the following equations, e must be read as € > e -*■ € 


EXAMPLE: 


[ A l e (1 ' p l e >E 3 e M l e j3 e ] ■ ( A l, a ' P le )e 3 e H le33 e ) + ( A l e (1 ‘ p l e )<; 3« H l,j3 e ) 
L J S,L S L 


D n *(m 1 )( a i c e i c {< 1 ' p i c ,h i c ji c " 1> + < 1 " p i e )H i e n.' 1) 

+ \ (M 1 )£ 1 S U1' 1 + ?2 9 e "ivVnajAju.l ) s, L 

0 

D 12 ■%)( A lc a ' P lc )e2l ' 1 c32 + A l e (1 ' P l e )e 2 H l e j2 ) S L 
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D 13 ' (i^) ( A l e <1 - p l e )e 3 e H l eJ 3 e ) 


S,L 


D 14 " (M°)( A l c (1 ' P l c )e 4 H l c jA ) 


S,L 


D l5 " (M^)( A l a (1_P l )e 5 T n G 


S,L 


D 16 ' (^) (*i a a ' P l )€ 6^ H ll»6 H lj nj 

S,L 


D 17 -° 


D 18 ' (57) ( A l a <1 - p l )U 8 b H lj8 b «8 t H lj8 t +? l 2 l 5H lj 11 a , ll b n) 1 b (e 8b ,i y8 b 4 *8 t I il aj 8 t 
D lll" (i^) ( A la < 1 - p l ){el i H 13 i ^u 5a iJ 11 a” r b» lt < i e: i bVu.'W 1 ) 

S,L 


n - 0 « n “ D 
D 1 10 u U 1 11 19 


D 21 " (1) ( A 2 e (1 " P 2 ) E l e H 2jl e ) 


S,L 


°22 “ (hJ ( A 2 e £ 2 { (1_P 2 ) H 2j2 " l} ) 


* 2 ' ' e 


S,L 
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D 23 " (m 2 )( A 2 e (1_P 2 ) e 3 e H 2j3 e ) 

S,L 

D 24 " °25 “ °26 " D 27 “ D 28 " °29 “ °2 10 " 0 “ °2 11 

D 31 “ (^2 )( A 3e (1_P 3e > % \& e ) 

J S,L 

D 32 " (Ml) ( A 3 e (1 ' P 3 e ) £ 2 H 3 e j2 + A 3 d (1 ‘ P 3 d ) £ 2 H 3 d J2 ) 

J 7 S,L 

°33 “(M^)( A 3e C 3 e {(1 " P 3 e )H 3 e j3e" l} + A 3 d *3 d { (1 ~ P 3 d )H 3 d j3<f l} 

+ V3 {(1 - P 3 )H 3J3- 1} ) 

S,L 

°34 " D 35 " °36 " ° 

D 37 " (M^)( A 3 d (1_P 3 d ) £ 7 H 3 d j7 + A 3 f (1-p 3 ) ’ C 7 H 3j7) 

S,L 

D 38 " °39 " D 3 10 = 0 " D 3 11 

D 41 “(w^)( A 4 c ( 1 “ P 4 ) £ l c H 4jl c ) 

S,L 

D 42 " (m^ ) ( A 4 c ^ 1_P 4^ £ 2 H 4 j 2 ) 


15 



0 


D 44 ■ ( 5 ;) (\ c C 4 t(1 “ p 4 ) H 4J« ' 11 '* W (1 ' P « ),, 4J4' 1, ) S>1 


°46 “ °47 ■ °48 ’ D 49 ' D 41<f D 4 11 ‘ ° ' D 45 


d 5i - (if) ( a s c i <i_p 5 ) yvvJ 


(^X A ‘ <1 ' P ‘ )C 5 H M5) s l 


D 52 ' D 53 ‘ ° ' ”54 ' (>f)( A 54 (1 - p 5 )i: 4 H 5J4) SJ 


A 5 £ 5 { <1_p 5 > H 5m5 _1 + ^1 (1 ' P 5 )G "ll a j“ a \™5 


+ A 5 E 5 {(l-p 5 )H 5j5 -l} 


' (if ) ( A 5 £ 6 { <1 " p 5 ) H 5m6 + Yl (1 ' P 5 )G ’W 


d 57 -° 


(if) (V“ (1 - p 5 )5 » 5 Ji b (E e b R Ual 8 b + 


D 511-(4) (V^S* Vs^i + %* H il a il a (e >l a + 
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D 5 10 ■ 0 * 

D 5 11 " °5 9 

D 61 "(m^)( 

i A 6 e l T ll (1-P 6 ) ^ ^fejl H 6w U b) 


S,L 

D 62 " D 63 " 

V ■ ° 

MO 1 

S E 5 (1 - P 6 ){H 6»5 + ? U E V? a I, 9 a J lI a \^) 

S, L 

°66 

[ A 6 e 6 { (1 " P 6 ) H 6m6 + T 11 (1 " P 6 )5 H 6nl^ Vjllg ^fernS -1 *) 

C 

D 67-° 


D 68 ‘((^) 

(Vll (1 ' p 6 )5 H 6 th1\j (E 8 b H iy8 b + % 


S,L 

D 6U‘ (i£ )' 

( A 6 (1 ‘ p 6 )H 6mll a {E ll b + T lf H lI a jll a (e iJa + H U b wll b )> ) 


D 6 10 * 0 * D 6 11 * D 6 9 


°72 “ (m^ )( A 7 d (1_P 7 ) E 2 H 7j2 ) 
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»»-(4)(V 1 ‘ P > ) E 3 d "7J3 d + A 7 t (1 - p 7> £ 

°74 * D 75 ' D 76 ’ 0 

D 77-(i^)( A 7 d E 7 d { < 1 -'>7> H 7 d J7 d - l} + A 7 £ *7 

°78 * °79 " °7 10 = ° = D 7 11 = °7 12 
D 81 " (m^)( A 8 £ 1 { (1_P 8 b )H 8 b jl + (1 " p 8 t )H 8 t jl 
+ W B l t Jl( (1 ' P 8 b )H 8 b Jl k 
D 82 ' D 83 ■ D 84 * ° 

D 85 ’( ) ( A 8 C 5 T H G "u b n>3 {<1 ' P 8b )H 8bjU a + ' 

D 86 ‘ ( ^ ) ( A 8 C i\° \ b J- a - P 8b >H 8bJll a + ° 


3 “733) 

S,L 

r {(1 " p 7 ) H 7j7 - 1} ) 

S,L 


+ a-VVl)*) 


1 ' p 8 t >H 8 t jn 0 1 ) 

S,L 

" P 8 t )H 8 t jna } ) 

S,L 
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D 88 ■(*:) (V (1 ' p 8 b He 8 b H 8bJ8 a + e 8 t \j8 t > + (1 ‘ p 8 t ) (E 8 b H 8 t j8 b + 


- e„ - e 


+ t 2 G H„ (e a . H . fi + e Q H , ft ){ (1-Po, )Hp 


8 b * e 8 t * T n° il b jll b V,t ’ 8 b ll a j 8 b T 8 t U a J® t y ' 8 b 8 bAia 


+ U'WjV 1 ) 


S,L 


D 811 " ( ST ) (V^ t(1 ‘ P 8 b )H 8 b #‘a + 


{1 + V ( \ b «ib + ^11 UW’) 


S,L 


D 8 10 " ° “ D 8 11 “ D 8 9 


“ill ■(% )( Vl H , te 31 {1 ' P11 a + ^ ( <1 - p U b ) \rfI b - 1 ) 

+ V»> - 1 ) 1 ) 


S,L 


\n “ D U3 " D i* 


D.., “ 0 


e 8 t H 8 t j8 t ) 
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Xs’ (ini) ( *11 5 H lb»5 { 1_P Ub + V ( a " P l)a ) "u a jU a _1 ) 


+ t 2 g h 


11 lfejlla 


f(l-p )H - l)}) 

V lib llb^b / ’ 


S,L 


■> 116 " (r) (‘SV^^Ub + V( a - p Ha )H u a m a '0 


+ ft *‘u a JU a (' 1 -Ollb ) %nS - l )>) 


S,L 


D U7 ‘ 0 


%8 ' (*[[) ( A u fc 8bV 8 > + E 8t H ua 78c ){l - P Ua + V^'W^b - *) 


+ T ? V J a ' P “a )H “A ' ') ’) 


S,L 


D n u' (>£) (*n { a ‘ P ua )H UaJll a ‘ lMe “. + ^ WS + T » 9 a 11 9 a j 9 a >) 

+ < d-fVX-Ub - lMe % + ^ ““.iVX + VnbViib*) 


S,L 


D u 10 " 0 “ D 11 11 
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D 10 1 " D 10 2 " D 10 3 “ D 10 A “ D 10 5 " D 10 6 “ D 10 7 " D 10 8 


D 10 9 “ D 10 10 " 0 " D io 11 " D 10 12 

D 11 1 “ D 11 2 “ D 11 3 “ D 11 A " D 11 5 " D 11 6 * D 11 7 " D 11 8 " 

D 11 9 " D 11 10 " D 11 11 " 0 " D H 12 


Variable 9: 

vf 0 

A<12 

m 2 

0 

Aq4 

ma 

0 

0 

Aq^ 

TT 7 
0 

0 

0 


+t q2, qA, and q7 are subscripts. 


Variable 10: 


BE ■ B s E s + Vl 
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<Vs,L-| 

f£) 

( A U a - p l )? T 10 5 H lj a 11 b mlo) 



S,L 

<V S ,L-' 

(£) 

( A 5 T 10 (1 " P 5^ H 5ml0 + T 9 G H 9 a j9 a H 9 b ml0 H 5m9 b *) 



S,L 

<Vs,L-< 

:v 

l ( A 6 (1 ~ P 6 )T 10^ H 6ml0 + T 9 G H 9 a j9 a H 9 b ml0 H 6m9 b * ) 



S,L 

^ b 8^S,L " ( 

Mg ) 

( A 8 T 9 T 10 G H 9 b ml0 { (1 ' P 8 b )H 8 b j9a + (1 " P 8 t )H 8 t j9 a ) ) 



S,L 

^ B 9^S,L “( 

,m 9 ) 

( A 9 T 10 H 9 b ml(/ 1_P 8 b + ? 9 G ( (1 - p 9 a )H 9 aj 9 a ' *) 


+ T 9 G H 9 a j9 a ( (1 " P 9b )H 9b m 9 b ' 


S,L 


All other elements are zero. 


Variable 11: 


b l ■ K l K 2 K 3 b 10 /M l 


10 

b 10 - (K l b 10 + n 10 >/M 10 * b 10 ■ (K? ’ V ' 
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All other bs are zero 


Without Filter Dome (total channel) 


Variable 12 


Dn " ( M i)( A lc e l C { (1_f i c )H l fc jl* - 1} + A l e e le {(1 ' P l e )H leJle 


f A la£l { d-P^ - 1 + Vl-; Pl )G H lj9a H v9b H 9fljl 


D 12 “(m^) ( A l c (1 " p l c )E 2 H l c j2 + A l e (1_P l e >E 2 H l e j2) 
D 13 " (M7)( A l e (1 ’ P l c )e 3 e H lej3e) 

D 14 “ (M^)( A l c (1-p l c )E 4 H l c j4) 

D 15 >B (M^)( A l a e 5 T 9 (1-p l )G H lj9 a H 9 b m5) 
D i6-(^)( A i a V9 (1 - p i )GH i j 9 a \ m e) 


D 17-° 



°18 “(m^) ( (1_p i )A l a {e 8 b H lj8 b + T U G H UUa (e 8t H Ub m8 t 

• + VSbVH. % aJ8b )} ) 


D 19 ‘ D 1 10 ' ° " ®1 11 ' D 1 12 
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(«:)( A 2 e a - p 2 )e i e H 2 3 iJ 


®22 * (m^ ) ( A 2 e t(1 ' P 2 )e 2 H 2j 2 ‘ e 2') 

D 23 '(ST )( A 2 e (1 ' p 2 )e 3 e H 2j3 e ) 


D„, “ D 25 = D 26 = D 27 = D 28 = D 2 10 = 0 ~ D 2 11 ‘ D 29 “ ° 2 12 


24 


D 31 ’ (*3 ) ( V 1 "**.* 1 . VO 

“32 - (b|)( A 3 e ( 1 - p 3 e )C 2 *^2 + H 3 d3 2 ) 

°33 ‘ (i ) - %' + (1 - p 3 d )H 3 d 13 d ' 11 

+ A3 £ e 3 {(l-P 3 )H, j3 - D) 
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°34 ’ 

°35 ‘ 

D 36“° 

°37 " 

( 4 )< 

, A 3 d (1 ' P 3 d )e 7 H 3 d j7 + A 3 f (1 " P 3 )e 7 H 3j7 

D _ m 

D_ 

= 0 ■ D__ ■» D_ ■ D 

i 

t U> 
00 

3 10 

39 3 11 u 3 12 

n m 

U 41 

(if) 



D 42 ' (m 4 ) (*4 c (1 - p 4 )C 2 H 4j2) 

°43 * ° 

D 44 ‘ ) (V‘ { W ' P 4 )H 4J4 - 11 + V4 ((1 - 5 4 )R 4J4 ' 1} ) 

D 45 ■ (fir) 

D 46 = °47 “ °4g " D 49 " °4 10 = 0 = °11 11 
D 51 -(h7)(Vi T 1 1 <1 - P 5 )G 

D 52 " °53 ’ 0 ’ D 54 " (mj)( A 5 (1 ” p 5 )£ 5 R 5j4) 

°55 ' A 5 £ 5 { (1 ‘ 5 >H 5m5 " 1 + T 11 (1 ' P 5 )G H 5mb "^“a 

+ V5 {(1 ' p i> fi 5J3 - 1>) 
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D 56 '(if ) (V6 (W S> H 5rf +, U° H 5-t *11.11. %*) 
d 57 -o 

D 58 ’ (if )( VW'V^t + T U G H 5mll b (e 8 b H ll a j8b + V’^ 

D 5 10 * 0 ' V * D 5 11 - D 5 12 
D 61 -(if)(VlV 1 - p 6 > <:H 6 ^ VlJ 

o 

»62 ■ D 63 • D 64 " 0 

D 65 ■(#) ( A 6 e 5 (1_P 6 ){H 6m5 +T 0° H 6«H Itfl \»5 } ) 

6' 

V 0 

D 68 " (if ) ( A 6 a ‘ P 6 )(£ 8 t H 6m8 t + V *6.1* ‘WdH + Vs^fejla V,*^ 1 ) 

D 6 10 ' 0 ’ D 69 ' D 6 11 ' D 6 12 
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72 


73 


74 


77 


78 


81 


82 


85 


(«7 ) ( A 7d W ~ P 7 )!: 2 H 7j2) 

(57 ) ( A 7d <1 “ P 7 )!! 3 <I H 7j3 d + A 7 f (1 " P 7 )£ 3 B 7j3) 

°75 - D 76 - 0 

(«7 )( V7 d { (1 - p 7 ,H 7 <i j7 d - 1} + V7 { (1 - p 7 )H 737 ' 1 
0 " °7 10 = °79 = °7 11 = °7 12 

(m^) (Vi a ' P 8b KH 8bJl + T U G W Vub W 

+ A 8 t £ l T u (1 ' P 8 t )G H 8 t Ml b ”u a 3l) 

D 83 ‘ D 84 " 0 

)( ViV^V 8 H 8 b 3ua V* + 

^ H 8 t m5 + T U G H 8 t mll b H Uajll a H U b m5^) 



°86 " (m 8 ) ( A 8 b e 6 Tll(1 ~ P 8 b )G H 8 b j u a ^m 6 + A 8t £ 6 (1 " P 8 t ) 

* H 8 t m6 + 1 il G H 8 tmllt} H ll a jUa 

D 87“° 

°88 = ( M^) ( A 8 b {£ 8 b ( (1 ' P 8 b )H 8 b j8 b ' 1 ) 4T ll (1 ‘ P 8 b )G XWVl^t 
+ V 8^1^115 **1^8^ }+A 8 t { e 8 t ( (1_P 8 t )H 8 t in8t" 1 ) 
■ K ll (1_P 8 t )G H 8 t n.ll i;> (e 8 b H n a j8 b + Si e: 8 t H llajna H llbm8 t ^) 

°8 10 " ° “ °89 = °8 11 °8 12 

D 10 1 = D 10 2 " D 10 3 “ D 10 A = D 10 5 = D 10 6 “ D 10 7 = °10 8 " 

D 10 10 " ° = D 10 9 “ °10 11 = °10 12 

D 11 1 = D 11 2 " D 11 3 " D n 4 = D 11 5 = °11 6 * °11 7 °ll 8 “ D H 9 

D 11 10 = D H 11 “ 0 “ D 11 12 * 0 
Bq = same as Variable 9 

Variable 13 ; 

B, •('£-}(*'. TT (l-p.)G 1. H . \ 

1 \ y \ la U l 12 1 ljUa Ub m 12/ 
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B 5 " (»Jr) ( A 5 T 10 (1 " P 5 ){H 5inl0 + T 9 G H 5m9 b H 9 a j9 a H 9 b mlO J ) 

B 6 “(m^ ) ( A 6 T 10 (1_P 6 )tH 6ml0 + T 9 G H 6m9 b H 9 a j9 a H 9 b mlO J ) 

B 8 “(m^) ( A 8 b T 9 T 10 (1-p 8 b )G H 9 b m9 b H 8 b j9 a + A 8 t T 10 (1_p 8 t ){H 8 t ml0 

+ T 9 G H 8 t m9 b H 9 a j9 a H 9 b ralO J ) 

All other elements are zero. 


Variable 14: 


b l " K l K 2 K 3 b 10 /M l 


b 10 ■ (K l b 10 + n 10 )/M 10 * b 10 ■ (l-H,) ’ 


10 


M 10 - 


All other b's are zero. 
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5.2 DYNAMIC SENSOR MODEL EQUATIONS 


Collecting all the above variables into a matrix differential 
equation, • 

Without Filter Dome: 


T - AT + DT + BpQ + BE + b + e 


( 8 ) 


v=C T T+C b b+n 


With Filter Dome: 


T = AT + DT + B.Q + B E +B T E T +b + e 

Q S S L L 


(9) 


v=C T T+C b b+n 


In the above equations (8) and (9), 


Eq. (8): 


T, T , Bq, B, b and e are the column vectors of 
9th order. 


Eq. (9):. 


A and D are 9x9 matrices. 

T, T , B~, B , B_ , b and e are the column vectors 
Q s L 

of 10th order. 

A and D are 10 x 10 matrices. 
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5 . 3 COUNT CONVERSION ALGORITHMS FOR M/WFOV TOTAL AND SHORTWAVE INSTRUMENTS 


A. Steady-State Linearization Algorithm 

1* Wide and Medium Field-of-View Total Instruments 

The algorithm contains two parts. First a table is generated off- 
line for each instrument. Using the dynamic model parameters computed, 
a table of gain and offset values is generated and stored. This table 
is updated as needed when new flight calibration data arrives. The second 
part operates on-line. As data from the instrument is recieved, it is 
processed by this portion of the algorithm to obtain the irradiance at the 
field-of-view. 

For the medium and wide FOV instrument the dynamic model has the form 


T * AT + DT + B q Q(T 2 ) + BE + b. 


v = C t T + C^b , 


( 10 ) 


4 th 

where the notation T designates the vector whose i component is 
T^; i.e., the fourth power of the i^ component of the vector T. The 
matrices and vectors A, D, B^, B and b are computed as specified in the 
model equations. 

The gain and offset values consist of elements of the form ( E Q ^» v o; ^» 
Y .), i = 0, 1, . . . , N, where 
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E , ■ iAE , i ■ 0, 1, . . ., N. 
oi o 


E q1 is computed by the above equation, while is obtained 
using the following algorithm using the function 


f (T ,E ) - AT +P.T BE + b 
o o o o o 


Algorithm for generating the table of points (Em , V, 


1. Initalize: E - iAE , T, - T . , z - 0(i-0 only), 

O O 1 ol-l 1 


k - 1, a - 1 


Compute* AT^, Vi 


4 ‘k ■ ^ b q <A + 40 *k >" 1 [f(T k> E) + e q *k’ 

«k - -a (A + AD I^r 1 [f(T k , E ol ) + B q (z k + A^)) 


Vl " z k + 


T, . , - T. + AT, 

k+1 k k 


The prime 


denotes the transpose. 
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3.-1 


a 22 is the (2,2) element of (A + 4D 
3 - If jh l f j (T k+l» E oi ) + ®Q 2 k+ll * j ' f j (T k’ E oi 


Set a - 2a/3 and GO TO 2 


«. if jij. lyVr V + b q Vi I 4 e , go to 7 


5. k - k + 1 


6. GO TO 2 


7 - V ol ’ C T Vl + c b b 

8. 6 » -CjfA + 4D t^)' 1 B 


V ■ Fqs Q 6 


9. Store E ol , » 0l> T ol 

10. If i > N, STOP 

11. i - i + 1 

12. GO TO 1. 

3 

The matrix is defined as follows 




where T^ is the ith component of T^. The algorithm generates the 
look-up table of gains and offsets off-line and stores them. 

The processing of data for count conversion is performed by the 
on-line part of the algorithm. Let V be a measurement to be processed 
consisting of readings of counts (v) , base temperature (T_) , field 
of view limiter temperature (Tp) and the copper heat sink temperature 
(Tjj) i i • 6. » 


V 



1. Obtain the value of v q ^ (table look-up) "closest" to v 
(i.e., |v - v q ^| is minimum). Note: start research at the last v q ^ 
used to minimize real time usage. 

/s 

2. Compute estimates E of irradiance at FOV 
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E - 


E oi + V (v - V oi> 

E ol + Y vi <v - v oi> + Y Bi (T B - T Bol ) 

+ Y F1 <T F ■ T Fol^ + Y Hi fT H ~ T Hoi* 


( 14 ) 


2. Wide and Medium FQV Short Wave Instruments 

The algorithm for the short wave instruments is similar to the 
total channel's, and consists of generating a table of gains and offsets, 
and an on-line algorithm for processing data. It is important to note 
that to obtain the estimate of the short wave irradiance at the field 
of view, the algorithm uses the estimate of the total irradiance corre- 
sponding to the some time to account for the effects of long wave 
radiation on the short wave measurement. 

The dynamic model of the medium and wide field of view short wave 
instruments has the form 

T = AT + DT 4 + B rt Q(T„) + B E + B- E- + b 
Q 2 s s 11 

V = C T T + C b b 

The table consists of elements (E sl , E^, V^, Y oim , 3 oim ) and 
uses the function 

f(T, E , E ) - AT + DT 4 + (B -B )E + B E + b. 0-5) 

s t S X S X t 
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Algorithm for Generating Table (Shortwave) 


1. E t ■ mAE t , z 1 - 0 (m = 0 only) 

2. E g - iAE g> T x - T ol _ r k - 1, a - 1 

3 - L \ ' 7 ^ V <a + 40 ? k rl tf( V V V + b q 
Vi * z k + Bz k 

4. AT k - -a (A + 4D i 3 )' 1 [£(I k , E s , E t > + B Q Vl ] 




5. 


6 . 


7. 

8 . 
9. 

10 . 

11 . 


Set 2nd element of T k+1 to T 2NQM ; i.e., T^ +1 2 ** t 2 nom 

u j*i ivw v V + b q vii 1 A |f j (T k- V V + b q z k 


Set a = 2a/ 3 and GO TO A 


I£ A I VW V E t> + B q Vi I 4 £ - G0 T0 10 


k - k + 1 


GO TO 3 


V oi« ■ °T Vl + c b B 
E oi» * - C T <A + 4D fy' 1 B 1 
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12 . 


8 - -C T ( A + AD T^) _1 (B s - B 1 ) 
Y oim " e'Q8 Q 8 


13. 

Store E , E , 
s’ t’ 

^oim* Y oim * 

14. 

If ii H GO 

TO 17 

15. 

1 = 1 + 1 


16. 

GO TO 2 


17. 

If m * M, STOP 

18. 

m = m + 1 


19. 

GO TO 1. 



The table is thus generated by the algorithm given above. This 
algorithm is used off-line to obtain the gains and offsets used in the 
on-line part of the count conversion algorithm. 

The processing of arriving "input" data from "flight" operations 
for count conversion may be performed by the on-line portion of the algorithm. 
Let V be a measurement of counts and base, FOV limiter and copper heat sink 
temperatures obtained from a short wave instrument; i. e.. 
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Arriving data of this form is processed as shown below. 

On-line Count Conversion Algorithm (Short Wave) 

1. Round E to the closest value of E . Note: E is the last 

tm 

output of the on-line count conversion algorithm for the appropriate 

A 

total channel. Let m be the index of the rounded value . 

tm 

2. Obtain the value of (table look-up) closest to v; i.e., 

|v - v . I is minimum for i = 0, i, . . . , N. Note: m is fixed by 1; 

1 oim 1 J 

start search at last value of v . used. 

oim 

A 

3. Compute estimate of irradiance at field of view aperture E g , 

A 

E v 

K - Ki + Ki* [ <V ' V oi m > ' 6 oln, (i - M, 

A A /A 

E = E - E 
1 s 

The on-line count conversion algorithm thus provide estimates of 
the short wave, long wave and total radiation incident on the field of 
view aperture in "normal" operation. The model parameters and look-up 
tables will be updated as needed when calibration data is obtained to 

account for changes in the instrument behavior. 
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Modified Kalman Filtering Algorithms 

1. Medium and Wide FOV Total Instruments 

The Steady State Linearization algorithms given are based on the 
assumption that temporal and spatial variations in the incoming radiance 
are sufficiently slow and small so that the output of the instrument 
is mainly determined by its steady state response. The algorithms 
presented in this section are based on concepts used in Kalman filtering 
and do not make this assumption so that the portion of the output due 
to the transient response is accounted for in the conversion. Further 
error reducing properties of these algorithms make them suitable for 
high accuracy measurements. 

The algorithm is recursive and consists of the following equations 
given in matrix form. 


v(t k ) 

- V(t k> - C T * (t k> - C b b 

(16a) 

i(t k > 

- f(t k ) + g t v(t k ) 

(16b) 

E(t k ) 

" *< T k-l> + «Vl> + °E V<t k> 

(16c) 

S(t k > 

- + C s V(t k> 

(16d) 


T(t k+1 ) = + ^ 4(t k ) + r Q Q<Vt k )) + r E(t k ) + 0S(t k ) + (ipb) 

where t k is the k c ^ sampling instant, At the sampling interval, V(t k ) 
the measurement vector at time t k consisting of counts, base, FOV and 


(16e) 
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heat sink temperature readings, and E(t k ) is the estimate of the 
irradiance incident on the FOV aperture at time t^. 


q, if T 2 (t R )< T 2N0M - e or ^T 2 (t fc )< T 2nom + e and Q ( T 2 ' 


o, otherwise 


The initial values of T(t Q ), E(t Q ) are obtained using the SSL algorithm 

A 

for generating the look-up table, and S(t Q ) = 0. 


2. Medium and Wide FOV Short Wave Instruments 
The algorithm for the short wave instruments uses the estimate 
E(t^) and S(t^) obtained from the algorithm for the total channels. 
This allows the algorithm to estimate effects such as the heating of 
the dome filter and leakage. 


v(t k ) = V(t k > - C T T(t k > - C b b 

(17a' 

T(t k ) • T(c k ) + G t v(t k ) 

(17b) 

£ s (t k > ■ Vw + ic K < w-i > + c e v<t k ) 

0 7c) 

S s (t k> ' S .<W + G S V(t k ) 

(17d) 

T(c k+1 ) - 4>T(t k ) + W 4 (t k ) + r Q P(f 2 (t k )) + r 8 E s (t k > 


+ + r T J(t k> + e T S <‘k> + *» 

(17e) 


A . 

E s^k^ is t * ie estimate of the short wave irradiance incident on 
the FOV aperture at time t^. 
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NON-SCANNEK DYNAMIC MODEL 


NODES 



H HEATER 

P TEMPERATURE PROBE 


Pif-nro 1 


ACR CONTROL ELECTRONICS 



Resistance 

Winding 






IU H- 
It03 
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Figure 







FUNCTIONAL FLOW DIAGRAM FOR SENSOR 
CHARACTERIZATION AND CALIBRATION 



Figure 4 
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Figure 







INFLIGHT MEASUREMENTS 




Figure 6 
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ORIGINAL PAGE IS 
OF POOTT ° T T M>TTY 

NONSCANNER THERMAL MODEL - WFOVSW 

(see List of Subscripts, p. vi, for explanation 
of nomenclature) 



8 APERTURE 

9 SW FILTER DOME 
H HEATER 

P TEMPERATURE PROBE 


Figure 7 
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MONSCANNER WFOVSW ENCLOSURES 

(see List of Subscripts, p. vi, for 
explantation of nomenclature) 


E 


FOV 
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sensor f ield-of-view (FOV) limiter into data counts, derives from this model two 
algorithms for the conversion of data counts to irradiance at the sensor FOV 
aperture and develops measurement models which account for a specific target source 
together with a sensor. The resulting algorithms are of the gain/offset and 
Kalman filter types. The gain/offset algorithm was chosen since it provided suf- 
ficient accuracy using simpler computations. 
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